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INTRODUCTION 

Linear carbon chains have been proposed as carriers of the diffuse interstellar bands 
(DIB) as long as 20 years ago (Douglas 1977). Recent experimental results by Fulara et. 
al. (1993) give new arguments for this interpretation of the origin of the DIBs. These 
authors show that several DIBs can be attributed to the species containing carbon and 
hydrogen atoms only. In particular, they suggest that the absorption features observed at 
6077 A and 7728 A correspond to the difFuse interstellar bands centered at 6089.7 and 
7721.7 A and attribute them to the 2 II- 2 n transitions in C 6 H and C S H, respectively. 

Recently Watson (1994) investigated the possibility that the linear carbon-containing 
chains may be the carriers of DIBs. Watson used a particle-in-the-box model to investigate 
the electronic transitions in homologous series of carbon-chain molecules. 

In the C271H chains, C 2 H and C 4 H have 2 E + ground states, with electron configurations 
•••(« — l) 7 r 4 n 7 r 4 ( 4 rc + l)cr. Two low lying excited states of 2 n symmetry are described 
by . . . (n — l) 7 r 4 n 7 r 3 ( 4 n + 1 )<j 2 and . . . (n — I) 7 r 3 7 i 7 r 4 ( 4 r 2 + l)cr 2 configurations. Longer 
chains have the 2 II ground state, with the configuration . . . (n - l)x 4 n 7 r 3 ( 4 n -f l)cr 2 
(Pauzat & Ellinger 1989; Pauzat et. al. 1991), so the (strong and parallel) 2 II- 2 n 
transition leads to the . . . (n — l) 7 r 3 n 7 r 4 ( 4 n + l)<r 2 configuration. We report here the 
results of an ab initio study of the electronic spectrum of C 4 H and relate it with the DIBs 
problem. 

METHODS 

The basis set used here is of valence tnple-zeta (Dunning 1971), double polarization 
(Dunning 1989) quality and consists of 110 contracted gaussian type orbitals. All ge- 
ometry optimizations were performed at the complete active space SCF (CASSCF) level, 
correlating one cf and eight 7r electrons on nine orbitals, using program GAMESS (Schmidt 
et. al. 1990). The geometry optimizations were restricted to the linear geometries. All 
four inter-atomic distances were optimized. 

The multireference configuration interaction (MR-CI) calculations were performed for 
several linear geometries of the molecule using MRD-CI method (Buenker & Peyerimhoff 
1974; Buenker 1980; Buenker 1982). In those calculations 17 valence electrons were 
correlated within 86 molecular orbitals. The standard selection/extrapolation procedure 
with Davidson correction (LanghofF Si Davidson 1974) was used. The resulting energies 
are denoted MR-CI(+Q). It has been shown (Knowles et. al. 1990) that this treatment 


85 


PMOIKMW 





Table 1. Optimized bond lengths of the X 2 S + , A 2 !!, and B 2 II states of C 4 H. 


State 

Method/basis 

Ci-C 2 

A 

c 2 -c 3 

A 

C3-C4 

A 

C,-H 

A 

B 

MHz 

411. 

ref. 

X 2 £+ 

Experiment 





4758.6567 

a 

X 2 £+ 

CASSCF/TZP2P 

1.2095 

1.3854 

1.1820 

1.0546 

4785.0 

b 

A 2 n 

CASSCF/TZP2P 

1.2613 

1.3505 

1.2139 

1.0538 

4706.3 

b 

B 2 n 

a ( r .. 1 

CASSCF/TZP2P 

i: - 1 1 aaa\ 

1.2544 

1.3534 

1.2716 

1.0526 


b 


° (Guelin et.ai. 1982) 
6 this work 


Table 2. Term energies for the electronic states of C 4 H 



X 2 £ + eq. 

A 2 U eq. 

B 2 n eq. 

-4> •- 

State 

MR-CI(+Q) 

(eV) 

MR-Ci(+Q) 

(eV) 

MR-CI(+Q) a 

(eV) 

Debye 


o.oo 6 

— 

— 

0.93 

A 2 n 

0.35 

0.01 

— 

5.41 

b 2 n 

3.70 

3.34 

3.24 

8.83 

C 2 £+ 

a - 11 un / 

4.99 

1 r\ \ 

4.61 

__ _ I _ 1 * 

— 

1.09 


a all MR-Cl(-fQ) energies relative to -152.3367709 E H 
6 Dipole moment 


results in less than 1 kcal/mole ( 0.06 eV) errors in the estimation of the full-Cl limit. 
Use of a basis set of triple zeta - double polarization quality should lead to similar errors. 
Thus, the relative energies here are accurate to better than 0.1 eV, which translates to 
130 A at 4000 A. 

AB INITIO RESULTS 

We have optimized bond distances in the 

three lowest-lying electronicstates of C 4 H: X 2 £ + , A 2 II, B 2 II. Table 1 shows the bond 
distances obtained in this study. 

The optimized geometries were then used in the multireference Cl calculations of the 
vertical and adiabatic excitation energies for four lowest lying states of C 4 H. The results are 
shown in Table 2. The boldface in the table indicates the adiabatic excitation energies. In 
addition we have calculated restricted bending potentials for the terminal C and H atoms, 
to confirm the linearity of the molecules in the ground and excited electronic states. The 
bending potentials in three lowest electronic states have minima at linear conformations 
of the molecule. 

The ground state of C 4 H is X 2 £+, but the A 2 II state is very low-lying, at only 0.01 
eV at MR-CI(+Q) level. The vertical excitation energy is 0.35 eV. The calculated energy 
difference of 0.01 eV does not give a conclusive assignment of the ground state. Kiefer 
et. al. (1992) do not obtain clear-cut assignment of the ground state of C 4 H and similar 
problems are discussed by Pauzat and Ellinger (1989). The rotational spectrum of C 4 H 
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Table 3. Theoretical band origins of C 4 H 


Transition 

Aoo (A) 

-a 1 uauu v/i igmd v/i 1 1 

transition moment (Debye) 

A-n-x 2 s+ 

infrared 

0.56 

B 2 n-X 2 £+ 

3884 

0.13 

C 2 £ + -X 2 £+ 

2484 

0.11 

B 2 n-A 2 n 

3869 

2.57 


has been observed in the 2 E + state, but this does not preclude non-zero population of 
the radiatively long lived 2 II state. On the other hand one may expect significant vibronic 
coupling between the E (A* in bent configurations) and II (A’ and A") states (Yamamoto 
et. al. 1987). This effect may contribute to the fact that our calculated B values for the 
2 E + and 2 II states bracket the results of astronomical observations. 

The vertical and adiabatic excitation energies between the B-X and B-A states differ 
by 0.46 eV ( 3.70 vs. 3.24 eV) and 0.09 eV ( 3.34 vs. 3.24 eV), consistent with the 
changes in the equlibrium geometries (see Table 1 ). 

Prediction of the band origins from ab initio requires estimation of the zero point 
energies of the involved states. Using harmonic vibration frequencies for the X and A 
states, reported by Kiefer et. al. (1992), the corresponding zero point energies are 5639 
and 5456 cm' 1 , respectively. In particular, the v z frequency (one of the triple C-C bond 
stretches) decreases from 2116 to 1864 cm" 1 , reflecting the weakening of this bond. 
Based on the increase in the C 3 -C 4 bond length, one may expect a similar effect on the r / 2 
frequency (other triple C-C bond stretch) in the B state, reducing its zero point energy. 
We will assume that the zero point energy of the B state will be reduced to 5250 cm -1 . 

The above information, together with the term values reported in Table 2. and the 
values of the dipole transition moment, allows for a prediction of the general features of 
the electronic spectrum of C 4 H. We will consider the following four electronic transitions: 
A-X, B-X, C-X, and B-A. The red end of the absorption spectrum will occur at the 0- 
0 transition energy and the maximum intensity should occur at the vertical excitation 
energy, corrected for the zero point energy of the lower state (provided it is larger than 
0-0 transition energy). The blue end of the spectrum is difficult to predict in the absence 
of detailed vibrational analysis, but one should expect relatively short progressions, as the 
changes in the bond lengths are not very large. The main parameters of the four electronic 
transitions are given in Table 3. 

Both B-X and B-A bands have a very similar origin A 0 o- A vibrational progression in 
the B-X system may occur because of the difference in the Ci-C 2 bond lengths of 0.06 A 
between the lower and upper states. In the B-A system, the Ci-C 2 bond lengths are very 
similar, but the C 3 -C 4 distances differ by 0.06 A. Consequently, a very similar vibrational 
progression may occur. However, the oscillator strengths differ by two to three orders of 
magnitude, parallel transition B-A being stronger. Thus, depending on the conditions and 
possible pathways to populate the very low lying A 2 II state, it may be possible to observe 
absorption from both states. Both bands should begin at around 3900 A and should 
overlap considerably. It is important to point out that the B-A band corresponds to the 
transition in the particle-in -the- box discussed by Watson. The C-X band, originating at 
2500 A may also be observable. 

PARTICLE-IN-TH E-BOX 

Chains longer than C4H have 2 II ground state (Pauzat & Ellinger 1989) and should 
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Figure 1 . Parallel bands in C 2 n H series and the DIBs 


have strong absorption spectra at low temperatures, between the . ..7t 4 7t 3 and . . . 7t 3 7t 4 
configurations. In a chain the ( n — l)x and mr molecular orbital have a nodal 

structure of the (n — l) th and n th eigenfunctions of the particle-in-the-box model. The 
electronic transition between those levels has been discussed recently by Watson (1994). 
In the case of a C 271 H chain, the equivalent box length and the wavelength of the transition 
between two lowest II states are given by the following relations: 

A(A) = 329.72 x l 2 {\)/(2n - 1) 
l{\) = ((2n- 1) x A(A)/329.72) 1 / 2 

where A is the transition wavelength and / is the effective length of the box. The wave- 
length of the parallel absorption band is changes linearly with the box length and conse- 
quently number of C atoms in a chain. 

C 4 H AND THE DIBs 

It has been proposed (Bell et. al. 1983), that the 4428 A diffuse interstellar band is 
carried by C 4 H. This assignment is consistent with a recent proposition of Fulara et.al.: 
Application of the simple model described above to extrapolate the proposed CgH and 
CgH lines to C 4 H gives the wavelength of around 4430 A. However, it is rather unlikely 
that C 4 H is responsible for the 4428 A DIB, as C 4 H is predicted here not to absorb above 
4000 A. Consequently, the 6077 and 7728 A lines seem to be too far into the red to 
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Table 4. Wavelengths of parallel transitions of C 2n H species, 


with n=2,3,4,5, and 6. 


Species 

box length (A) 

wavelength (A) 

C 4 H 

5.93 

3869“ 

c 6 h 

8.53 

4798 

c 8 h 

11.13 

5834 

CioH 

13.73 

6906 

Cj 2 H 

16.33 

7993 


a from MRD-CI calculations 


arise from (relatively) short CgH and C 3 H molecules and to be too far apart to arise from 
molecules differing by only two carbon atoms. 

However, it is interesting to predict the wavelengths of the transitions in C 2n H series 
by using the partide-in-the-box model. Our best estimate of the energy of the B 2 II-A 2 n 
transition in C 4 H corresponds to the wavelength of 3869 A. The equivalent box length is 
5.93 A. We will now assume that the effective box length for the C 2n H molecules is equal 
to the length of the box of the previous member of the series, augmented by 2.6 A. 2.6 
A is a reasonable estimate of the sum of the lengths of two adjacent C-C bonds in a 2 II 
state of the molecule ( see Table 1). 

Table 4. shows the predicted transition energies for C 4 H - Ci 2 H, using the box lengths 
described above. Figure 1. shows a plot of the diffuse interstellar bands observed against 
HD183143, with the data taken from (Jenniskens 1992). In the figure, the height of the 
bar is proportional to the equivalent width per color excess and the values from Table 4. 
are shown with thick lines, with an arbitrary intensity. 

The following coincidences can be observed between the predicted frequencies for 
C 2n H species and strong diffuse interstellar lines: (4798; 4761.67), (5834; 5779.48), 
(6906; 6939.00), and (7993; 7927.80). Obviously, the ab initio calculations can not 
achieve spectroscopic accuracy. Thus, the assignments shown above are quite tentative 
and serve more as an illustration than the genuine assignment of the lines, particularly 
considering the simplicity of the model used. 

Nevertheless, the RMS deviation of the assignment given above is only 50 A and 
is smaller than errors in ab initio calculations and comparable with the matrix shifts in 
absorption spectra. The coincidence of the calculated transition wavelengths with four 
of the intense diffuse interstellar bands is generally consistent with the idea that DIBs 
originate not from one, but from a small group of molecules (Watson 1994). Also, this 
interpretation of the origin of the DIBs is consistent with the lack of observation of DIBs 
in ultraviolet. 

A very similar reasoning should also apply to the C 2n „iN series, as the CH group is 
isoelectronic with the N atom. An ab initio study of C 3 N (Sadlej & Roos 1991) predicts 

the wavelength of the 2 II- 2 n transition at around 3600 A, consistent with a shorter length 
for the terminal C-N bond in C 3 N, compared to the C-C bond in C 4 H 
CONCLUSIONS 

We have presented an ab initio study of the low-lying electronic states of C 4 H. This 
species has a X 2 £ + ground state and two low lying II states. Based on the vertical and 
adiabatic excitation energies between those states we suggest that the 4428 A diffuse 
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interstellar band is not carried by C4H. The application of the particle in a box model 
shows strong coincidences between the strong DIBs and predicted wavelengths of 11-11 
transitions in C 27 M series. Based on those coincidences, we suggest the Q. 2 rM species as 
good candidates for carriers of diffuse interstellar bands. 
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